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ABSTRACT 



An analysis and test rig measurements of the per- 
formance of a transonic axial turbine are reported. The 
purpose was to confirm the accuracy of measurements made 
in a test rig which was designed to separate the losses 
occurring in the stator from the losses occurring in the 
rotor blade rows. The analysis was programmed for the 
Hewlett-Packard 21-MX computer. Reasonable agreement 
between predicted and measured characteristics was obtained 
using experimentally determined losses in the computer 
prograun. Lack of agreement was noted using theoretical 
values. It was concluded that the rotor was not choked 
at the conditions in the tests, and that the test rig 
measurements were valid. A successful technique for 
smoothing the data obtained from the rig is also reported. 
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I. INTRODUCTION 



The trainsonic tiirhine test rig installation at the 
Naval Postgraduate School Turbopropulsion Laboratory was 
designed to study the effects on turbine performance of 
varying axial and tip clearances, to study the effects of 
blading design on turbine performance, and to allow the 
separate determination of stator and rotor losses in an 
operating machine . 

Until the work of Solms (Ref. 1) the separation of 
rotor and stator losses through test rig measurements had 
not been attained satisfactorily. Through improvements 
in hardware and instrumentation and improvements in the data 
reduction process the stator and rotor losses were deter- 
mined separately and were reported in Ref. 1. 

Anomalies remained, however. Specifically, the turbine 
coitiLfdguration designated as Turbine C in Ref. 1 gave 
different results when compared in terms of "referred" 
quantities depending on whether the discharge was to 
atmospheric pressure or to a region of reduced pressure. 
(Turbine C had converging-diverging stator passages in an 
axial entry, single impulse stage. The turbine was designed 
to operate in the transonic range.) In addition, R&f. 1 
reported considerable scatter in the loss coefficients. 

In the work of Robbins (Ref. 2) it was shown that 
discharge pressure affected the measurements of flow rate 
into the stage and also affected the labyrinth leak rate. 
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Accordingly, Robbins determined accurately the flow rate 
into the stage and the leaJc rate through the labyrinth seal 
(See Fig. 1) for all operating conditions. . The results 
then obtained for turbine C were reported fully in Ref. 2. 
The continued presence of scatter in the measured loss 
coefficients was reported and a smoothing technique to 
eliminate the scatter was suggested. 

Before the Turbine Test Rig could be used to measure 
the effects of varying parameters i 

1. The scatter in the loss coefficients had to be 
eliminated. 

2. The overall accuracy of the performance results 
evaluated from the rig measurements had to be 
verified in some way. 

The resolution of these problems was the goal of the 
present work and is the subject of this report. First, 
a satisfactory method was found for smoothing the loss 
coefficients. The method is described in Section III. 

The approach taken to verify the performance of the 
rig Was to first devise an analysis which predicted the 
performance of the turbine in terms of unknown loss coeffic' 
ients, and then to show that the measured loss coefficients 
were consistent with the predicted behaviour. 

A description of the Turbine Test Rig is given in 
Section II. The analysis of the behaviour of the turbine, 
involving a comparison of an analytical prediction with 
the results of a short test program, is described in 
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Section IV. Details of the analysis and the computer 
program are given in Appendix A. 
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II. TURBINE TEST RIG INSTALIATION 



A. DESCRIPTION 

The test installation consists of three major compon- 
ents; an Allis-Chalmers twelve stage axial flow compressor, 
an exhauster assembly, and the turbine test rig (TTR) 
itself. 

The compressor is the source of driving air for the 
TTR and for the exhauster assembly. Fig. 2 shows the 
piping arrangement. Turbine air passes through the first 
settling tank into an eight-inch pipe containing a flow 
nozzle, into the second settling tank and into the turbine. 

Fig. 3 shows the plenum, the floating stator assembly, 
the rotor, and the dynamometer (Ref. 1). Pressure ratios 
of 6*1 can be achieved when the system is hooded. The 
hood was needed to achieve high pressure ratios in the 
tests reported here. Fig. 4 shows the turbine blading 
of the stator and rotor. Ref. 3 contains detailed descrip- 
tions of the test rig hardware. 

The floating stator assembly shown in Fig. 3 permits 
measurements of the axial force and the torque on the 
assembly. Axial and rotational movements axe constrained 
by calibrated force transducers that axe heat insensitive. 
These measurements, together with wall static pressure 
measurements, allow the determination of the average 
axial and tangential velocity components at the stator exit. 

In this report one configuration designated Turbine C 
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Was tested, the geometry of which is shown in- Fig. 5* 

Table I describes the geometry quantitatively. The stator 
blade profile is shown in Fig. 6. The blades of the stator 
generate a converging-diverging nozzle shape. Pressure 
measurements were taken at the locations shown in Fig. 6. 
The pressures necessary to the analysis of the stator axial 
force were taken at the locations shown in Fig. ?• 

B. TEST MEASUREMENTS AND ACCURACY 

1 . Mass Flow Rates 

Appendix A of Ref. 2 gives a detailed description 
of the method used to determine both the turbine flow rate 
and the labyrinth seal leak rate. 

2. Forces. Torques. Temperatures, and Pressure 
Ref. 3 and Ref. 5 give calibration procedures 

for the TTR. Identical procedures were employed here. 

Table II of Ref. 1 gives the expected accuracies of the 
me asurements . 

C. TESTING AND DATA REDUCTION 

The TTR data collection system is described in Ref. 4. 
Appendix D of Ref. 1 gives a detailed explanation of the 
turbine test procedures. Those procediares were followed 
here with the exception that a constant RPM was held and 
the pressure ratio varied over the desired range. The 
data reduction method developed in Ref. 1 and Ref. 2 was 
revised as described in Appendix B. 
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III. DATA SMOOTHING TECHNIQUE 



The sensitivity of the loss coefficients to variations 
in measxired quantities was shown in Ref, 2. It was also 
stated that some variation is unavoidable since measure- 
ments are taken over a period of more than a minute. It 
was also pointed out in Ref, 2 that the parameter most 
important to the calculation of the loss coefficients was 
(the average pressure at the stator exit). is an 

average pressure that can not be measured directly. It 
is derived as described in Ref. 1 from many other measure- 
ments, Significant scatter was observed in the variation 
of P^ as speed was varied at fixed pressure ratio. However, 
it was found that the hub and tip pressures (P^^ and P^) 
measured just downstream of the stator varied smoothly over 
the same range. Since these two pressure were measiired 
directly, and since they varied smoothly, it was assumed 
that the pressure behind the stator must vary smoothly also. 
Consequently, it was determined that a polynomial curve fit 
could be used to describe the variation of P^ as speed was 
varied. The variation of P^ was represented as a function 
of P^, P^, and the isentropic head coefficient (K^g) 
in the formi 

^ — . ( 1 ) 

Where is the polynomial coefficient. 

Then, (where P^^ is total pressure upstream of 

the stator) was computed using the expression 
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Fig. 8 is an illustration of this procedure, which 
was added to the "bulk process" data reduction program. 

The data points in Fig. 8 were taken from Runs 6 and 7 
in Ref. 2. The two Runs were for the ssane conditions but 
were made at different times. It can be seen that the 
trends are the same for both runs but the scatter is 
considerable. The scatter is the cause of the scatter in the 
calculated loss coefficients. The lines on Fig. 8 are the 
polynomial approximations according to Eq. (1) for the two 
runs. It can be seen that the polynomial approximation 
averages the data and maintains the original trend. The 

S "t 

chosen smoothing function was a 1 degree polynomial. 

Fig. 9 shows the stator loss coefficient prior to 
smoothing for the points in Run #?. Fig. 10 shows the 
same data after smoothing. 

Similarly, Fig. 11 shows the rotor loss coefficients 
for Run #7 prior to smoothing and Fig. 12 shows the same 
data after smoothing. 

It can be seen from these figures that the scatter 
has been removed. This is particularly true for the totor 
loss coefficients. As pointed out above, it is believed 
that the observed scatter was due to the sensitivity of 
the loss coefficients to small changes in measured quan- 
tities during the data collection process. The smoothing 
technique removes the random variations recorded during 



15 



data collection and results in a much more realistic 
representation of the variation in the losses. 

The smoothing technique was incorporated into the 
data reduction program as described in Appendix B. 
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IV. ANALYSIS OF TURBINE PERFORMANCE 



A. APPROACH 

In order to determine if the performance evaluated 
from the rig measurements was accurate an analysis to 
predict the behaviour of the test turbine was carried out 
and programmed in BASIC language. 

A performance test was then conducted in a particular 
Way in order to provide a comparison of the measured with 
the predicted behaviour. 

I 

B. ANALYTICAL PREDICTIONS 

The flow through the turbine was analysed using a 
pseudo-1 Dimensional compressible approach. The analysis 
is described in detail in Appendix A, together with the 
computer program which was used to obtain predictions of 
the turbine performance. One of the inputs which the 
program requires is the rotor passage loss coefficient at 
zero incidence. Using the method given by Vavra in Ref. 6 
it was determined that the rotor loss coefficient should 
have the value .251^. However, results of previous tests 
of the turbine indicated that the rotor loss coefficient 
was rarely as high as . 251 ^ and could be as low as .1. 
Therefore, the performance of the t\irbine was analysed 
using rotor loss coefficients of .2514 and .1. The pre- 
diction program was run for both 15,000 and 18,000 RPM with 
an assumed rotor loss coefficient of .2514, and for 18,000 
RPM with an assumed rotor loss coefficient of .1 
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The above parameters were chosen to obtain a prediction 
of the turbine performance in a range in which experimental 
data could be obtained. In particular, the analysis could 
be used to predict the pressure ratios at which choking 
would occur in the rotor as well as in the stator. The 
pressure ratio at which choking occurred in the rotor could 
then be established experimentally in a test conducted at 
fixed speed. An examination of the choking condition was 
considered to be a first test of the performance analysis. 

The results of the analysis for the parmameters given 
above are shown in Fig. 13. Fig. 1^, and Fig. 15. What is 
shown in the figures is the map of the range of values 
which unknown parameters in the analysis can have, that 
leads to a solution. 

Figures 16-19 show predicted performance parameters 
for the case of 18,000 RPM and assumed rotor loss of .1. 
These results will be discussed in conjunction with the 
results of the turbine test run. 

C. EXPERIMENTAL TURBINE TEST 

In order to examine the occurrence of rotor choking 
the turbine was run at constant RPM and the pressure ratio 
across the stage was increased in increments by lowering 
the back pressure. The point at which the horsepower ceased 
to increase for an increase in pressure ratio was examined 
to determine the choking point. At the condition where 
the flow reaches a Mach number of unity at the exit of the 
rotor, the power produced by the turbine can not be changed 
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by altering the downstream pressure. 

Tests were conducted in this manner at 15 t 000 and 
18,000 RPM. The controlled parameters for the tests are 
given in Table II. The reduced data is given in Table III. 
The referred horsepower is shown plotted versus the pressure 
ratio for the 15,000 RPM run in Fig. 20 and for the 18,000 
RPM run in Fig. 21. 

The results are discussed in the next section. 
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D. COMPARISON OF ACTUAL AND PREDICTED PERFORMANCE 

Figures 13 f 1^» and 15 are maps of possible solutions 
for the flow through the turbine, when the parameters which 
are unknown are allowed to vary. It is noted first that 
the predicted range of solutions extends to pressure 
ratios below the predicted choking line. The 

explanation for this apparent anomaly is that the program 
calculates all possible solutions, and for any point below 
the line a particular combination of losses and blockage 
factors existed which would allow a solution at that point 
without choking. It should be pointed out that any point on 
the plot can be brought to the choking' line by reducing 
the stator loss by a very small amount. This is the 
procedure that was followed to get the range of values 
at choking shown in Figures 16-19 . 

Fig. 13 is the map produced by the program for RPM 
equal 15 » 000 and for an assumed rotor passage loss coeffic- 
ient equal to .251^. As can be seen from the figure the 
pressure ratio, choking was about .28, cor- 
responding to a stage pressure ratio 3* 57* 

Fig. 14 is the map for RPM=18,000. Note that the predicted 
stage pressure ratio for choking is aigain about 3«57» 

Fig. 20 and Fig. 21 show the variation of the referred 
horsepower vs. the pressure ratio which were measured in 
turbine tests at 15 t 000 and 18,000 RPM respectively. It 
can be seen that the referred horsepower did not become 
independent of the back pressure at the pressure ratios 
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achieved in these tests. It was concluded that the stage 
did not choke at the predicted pressure ratio of 3 •57* 

It was noted however that at the highest pressure ratios 
obtained at 18,000 RPM, the slope of the horsepower curve 
was becoming smaller. 

The velocity diagrams in Figures 23-31 are also consis- 
tent with the argument that the rotor did not choke during 
the turbine tests. The velocity from the rotor, 
behaves smoothly and in a predictable manner throughout 
the pressure range at both test speeds. This might not 
be expected if, following choking, shock waves appeared 
downstream of the rotor exit plane. 

Note the values of rotor loss coefficient given for 
the test results in Table III. With the exception of 
point #10, which was considered to be in error, they were 
all smaller than the value (.251^) assumed in the first 
performance calculations. As explained in Appendix A, 
the value of the rotor passage loss coefficient used in 
the program is the smallest value that can be calculated 
for the overall rotor loss coefficient. Therefore, if 
the experimental results were accurate, the computer program 
could not predict the correct performance since it could 
never calculate a rotor loss less than the input value of 
the passage loss, which was .2514. It was as a consequence 
of this observation that the performance was re-calculated 
using a rotor passage loss coefficient equal to .1. 

Figures 16-19 show the choking behaviour of the 
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turbine predicted using the computer program for an assumed 
rotor passage loss of .1. In these figures, the data for 
test point #6 is also shown. In this case it can be seen 
that the predicted pressure ratio (^2^^to^ choking is 
about .245, corresponding to a stage presstire ratio (P-|;q/P 2 ^ 
of about 4.08. The latter is slightly higher than the 
maximum pressure ratio that was attained in the turbine 
tests. (At 18,000 RPM the highest attainable pressure 
ratio was 3.97 t due to the characteristics of the dynamom- 
eter.) 

Fig. 17 is the range of efficiencies predicted for the 
choked condition. It is noted that the efficiency measured 
at test point #6 was reasonably close to the predicted 
maximum efficiency. 

Fig. 18 shows the range of stator losses for which 
solutions existed, in comparison with the value of stator 
loss measured at test point #6. It can be seen that the 
measured value intersects the predicted range of possible 
solutions. Fig. 19 shows the predicted range of the rotor 
coefficient. Again, it can be seen that the measured value 
at test point #6 overlaps the predicted range of possible 
solutions. 

Fig. l6 shows the predicted range of horsepower 
compared with the horsepower measured at test point #6. 

The measured horsepower was slightly greater than the max- 
imum value which was predicted. 
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E. DISCUSSION 

The results shown in Figures 16-19 illustrate the 
uncertainty in the prediction of the performance of the 
turbine when only the rotor passage loss and blockage 
factor are known. It is recalled that the analysis 
satisfies only continuity through the stage, and values of 
an additional loss coefficient and an additional blockage 
factor must be established to obtain a unique solution. 

It is the aerodynamic shaping of the sxirfaces which deter- 
mines these factors. Figures I6-19 show the possible 
range of performance that results simply from the areas 
of the passages and the blade angles. 

The irregular shapes of the bounds on the possible 
solutions in Figures I6-19 are interesting. There is no 
obvious explanation for the reduced range of solutions near 

kbi = .81. 

The results obtained in the present work suggest 
that the performance of the turbine as measured in the 
turbine test rig is reliable. It had been thought prev- 
iously that the rotor loss measurements were too low. 

Here, an analysis was carried out to predict the performance 
of the turbine, and the predicted results have shown very 
good agreement with the performance measurements in the test 
rig. 

Whether or not the rotor was choking had also been a 
question in the past. In the present work both the analysis 
and experimental tests have shown that the rotor was not 
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choked in any test to date. The computer program with an 
assumed rotor loss of .1 has shown that the test rig should 
not be choked at any pressure ratio below about 4.08, Test 
data has shown that the rotor was not choked at a pressure 
ratio just slightly below 4.0. 

It is also of interest to note that for a rotor 
passage loss coefficient of ,2514 the predicted pressure 
ratio was shown to be independent of RPM, At both 15 #000 
and 18,000 RPM the choking pressure ratio 
calculated to be about 2 )> 57 ^ More results are needed to 
confirm that the choking pressure ratio is indeed indepen- 
dent of speed. 
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V. COxNCLUSIONS AND RECOMMENDATIONS 



1. The rotor of Turbine C was not choked at any pressure 
ratio tested so far. This conclusion is based on the 
prediction of the computer program and also on the 
turbine test results. 

2. The loss coefficients measured in the TTR and smoothed 
as described in this report can be accepted as being 
truly representative of the losses in the separate 
blade rows of the turbine. 

3 . In particular, the results from the computer program 
suggest that the magnitude of the measured rotor loss 
coefficients is probably correct. The method used 

by Vavra in Ref. 6 predicts much larger values of loss 
coefficients for the present rotor geometry than those 
which were measured. When the measured value of the 
rotor loss coefficient was entered into the computer 
• program, there was good agreement between all the 
calculated and the measured turbine performance 
parameters. When the higher (calculated) value of the 
rotor loss was entered, there was a pronounced dis- 
agreement between the calculated and the measured 
performance . 

4-. The computer program should be used (and progressively 
developed) in conjunction with all tests carried out 
in the turbine test rig. Most importantly, a test 
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should be conducted to determine the rotor choking 
condition experimentally. This will require the purchase 
of a water-brake dynamometer to extend the power-speed 
range of the test rig. 

5. Since the accuracy of the test rig measurements is 

no longer in doubt, experiments to determine the effect 
of parameter changes (e.g. axial and tip clearances) 
on turbine performance can go ahead. 
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TABLE I 



TURBIHE GEOMETET 





TURBINE 


STATOR 


ROTOR 






C 


1 


1 






DESCRIPTION 


2 

(in? 


NUMBER 

OP 

BLADES 


“'(no 


h/ . X 

(in.) 




STATOR 1 


.CONVERGING 

DIVERGING 

NOZZLES 


2.9050 


31 


4.184 


0.5775 




ROTOR 1 


CIRCULAR 

ARC 


7.119 


60 


E„-4.19? 

H^j-4.250 


h^=0.732 

h2»0.8475 
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TABLE II 



PARAr/IETERS FOR TURBINE TEST 



POINT # 


RPM 


Pto/^2 


1 


15,000 


2.01 


2 


15,000 


2.^2 


3 


15,000 


2.93 


4 


15,000 


3.^9 


5 


15,000 


3.68 


6 


18,000 


3.94 


7 


18,000 


3.47 


8 


18,000 


2.98 


9 


18,000 


2.52 


10 


18,000 


2.02 
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POINT #1 



TABLE III TURBINE TEST RESULTS 



PRESSURE RATIO 


= 2.011859838 


REF FLOW RATE 


-• 1 , 020000000 


REF RPrI 


= 14262.15617 


REF ROTOR HOMENT 


- 102.5192926 


REF HP 


= 23.19933152 


EFF. T-S • 


= 0.712728663 


STATOR LOSS THEOR, 


= 0.107100795 


S TAT 0 R L 0 S S C 0 E F" f ■' . 


0.220295023 


ROTOR LOSS COEFF. 


= 0.198187353 


ROTOR LOSS THEOR. 


= 0.386101836 


ISEN. HEAD COEFF, 


= 4,1036270755 


T H . D E G . 0 F R E A C T 1 0 N = 0 „ 0 2 5 3 8 6 6 4 4 
A C T . DEG. 0 F R E A C T 1 0 N = - 0 , 0 8 4 6 7 8 1 8 2 


VELOCITY TRIANGLE 


DATA 


VI = 973. 73327 


V2 = 204,23594 


VA1=-- 226 1.65467 


VA2= 155.42281 


'ii'iji = 946, 98688 


VU2= 132,49932 


ALPHA l--=' 76.5398 


:7 ALPHA 2= ^i0. 44789 



TH fel-J . 444 1 U B E I R 2 = ~ 5 9 84574 



POINT #2 

TresoIjf^e Rfrrid'“ 
REF FLOW RRTE 
REF RPM 

REF ROTOR nOMEHT 
REF HP 
EFF. T-3 



STATOR 
STATOR 
ROTOR 
RU FOR 



LO 

LO 

LOS 

LOS 



3S THEOR, 
3S COEFF, 
3 COEFF. 

5 THEOR. 



•-I .< ••j |-j c '7 1-1 • -1 

ill . -t ill O ■_! i O O O 
1 , 020000000 
14309. 05 S8 
1 3 2 . 8 8 4 / 9 F 
3011 04483385 
0. 748970950 

0. '1071 00795 
0, 183632754 
0. 133537254 
0. 455631440 



ISEH. HEAD COEFF. = 4.941765312 
TH. DEG. OF REACTION ==-8 . 76562E-03 
ACT. DEG 11 UF REACT I UN=-£i . 272 F9690 1 



VELOC 


ITY TRIANGLE 


D fl T f i 




VI = 


1 1 2 fo . U 6 8 4 




184.36934 


VA 1 •--= 


265 , 46949 


VA2= 


178.98521 


VU 1 - 


1094. 32844 


VU2 = 


4 4 , 2 3 01 6 6 


ALPHA 


1= 76.364 




ALPHA 2= 13.88078 


BETA 


1 = b 4 . W 2 -8 


13 


BETA 2 = -70.79960 



29 



TABLE III TURBINE TEST RESULTS (CONT.) 



POINT #3 



PRESSURE RfiTIO 
REF FLQW RRTE 
REF RPM 

REF ROTOR MOMENT 
REF HP 
EFF„ T-S 



= 2.932Sj?8?4 
== 1.020086000 
■-= 14297,06424 
= 161.7256587 
“ 36.6868547’ 

= 8 . F 6 5 2 ii b 6 4 



STATOR LOSS THE OR. 
S T Fi T 0 R L U S S L U E. F F . 
ROTOR LOSS COEFF. 
ROTOR LOSS THE OR, 



'0 . 1 07 1 00795 
II 1 r’ 2 0 1' 5 I’’ 8 r* 
0 . 093520223 
0,43179246 



I 'E'EN I. HEAD COEr'F. 5.882942922 

T H . DEG. 0 F R E fi C T 1 0 N = - 0 , 0 3 1 2 3 9 1 3 9 
A C T . D E G , 0 F R E A C T 1 0 fi -■= - 0 ,2417 9 5 5 9 2 



VELOCITY TRIANGLE DATA 



VI = 1248.51552 V2 = 
VA1= 314.71895 VA2= 
VU1= 1208, 19324 VU2= 



2 1 b . 9 8 7 6 3 
2016. 88375 
-71 .64356 



ALPHA 1- 75.39969 ALPHA 2= -19.10090 

BETA 1 = 64.49981 BETA 2 = -PI. 78607 



POINT #4 



PRESSURE RATIO 
REF FLOW RATE 
REF RPM 

REF ROTOR MOMENT 
REF HP 
EFF, T-S 



= 3.490951547 
= 1 . 020000000 
= 14290.27554 
= 181,5713712 
= 41.16922187 
= 0,753539626 



STATOR LOSS THEOR 
STATOR LOSS COEFF 
ROTOR LOSS COEFF. 
ROTOR LOSS THEOR. 



0. 107100795 
0. 158068259 
0. 134491135 
0.427987152 



I SEN. HEAD COEFF. = 6-703669585 
TH. DEG. ’OF REACTION =--0.014684530 
A C T . D E G . 0 F REA C T 1 0 N = - 0 . 2 1 3 8 8 4 7 3 8 

VELOCITY TRIANGLE DATA 



VI = 1333,15077 
VA1= 355.00633 
VU1= 1285 II 01 420 



V2 = 282,31378 
VA2= 238.26055 
VU2= -151.43638 



I 



ALPHA 1= 74.55631 ALPHA 2= -32.43969 

BETA 1 = 64.26912 BETA 2 = -71.41323 
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TABLE III TURBINE TEST RESULTS (CONT.) 



POINT #5 



PRESSURE RRTIO 
REF FLOW RhTE 
REF RPM ■ 

REF ROTOR MOMENT 
REF HP 
EFF. T--S. 

STRTOR LOSS THEOR. 
STATOR LOSS COEFF. 
ROTOR LOSS COEFF. 
ROTOR LOSS THEOR. 



- 3.677359655 
=- i. 020000000 
= 14292.77432 
= Is 6. 8682436 
= 42.3776356 
•= 0 . 7528013647 

= 0.107180795 
= 0. 140979499 

- 0. 192S6i’'968 

- 0.435423887 



ISEH. HEAD COEFF. 


:r: 


6.9505849: 


TH. :DEG. of REACTION = 


3. 90802E-I 


A C T . DEG. 0 F REA C 


TION= 


-0 , 2236967 


VELOCITY TRIANGLE 


DATA 




VI - 1362.18537 


V2 = 


301 . 4299 


VA1= 361.10549 


VA2= 


250 . 0096 


VUl= 1313.45034 


VU2= 


-168,3900 


ALPHA 1= 74.627 


rr -1 


ALPHA 2=-- 


BETA 1 = 64.690 


O •!' 


BETA 2 = 



~ 3 3 . 9 6 1 6 2 
--7 1.02 123 



POINT #6 



PRESSURE RATIO 


= 3.938522427 




REF FLOH RATE 


1 . 020000000 




REF RPM 


= 17091.82469 


1 


REF ROTOR MOMENT 


= 169.2232867 


i 


REF HP 


= 45.89160403 




EFF. T-S 


= 0.784488861 




STRTOR LOSS THEOR. 


- 0.107100795 


1 


STATOR LOSS COEFF. 


0.141347589 


1 


ROTOR LOSS COEFF. 


— 0 . 1 S i' 1'' 3 3 f' 4 0 




ROTOR LOSS THEOR. 


= 0.421689292 




ISEN. HERD COEFF. 


= 5.050902877 




TH. DEG. OF REACTION = 0. 184266156 
ACT. DEG. OF REACT ION= 0.175212937 




VELOCITY TRIANGLE 


DATA 


♦ 


VI = 1318.03669 


V2 = 269.74423 




VR1== 334.28819 


VR2^= 260,54759 




VU1= 1274.94005 


VU2= -69.83484 




ALPHA 1= 75.3078 


2 ALPHA 2= -15. 9043c 





BETA 1 = 61.53822 BETH 2 = -70.56818 



31 



TABLE III TURBINE TEST RESULTS (CONT. ) 



POINT #7 



PRESSURE RflTIG 


= 3.467558140 


REF FLOW RATE 


= 1.920000000 


REF RPM 


= 17039.50815 


REF ROTOR MOMENT 


= 155.9660642 


REF HP 


= 42.16691707 


EFF. T-S 


= 0 B 7 8 4 1 7 iS 4 4 2 


STATOR LOSS THEOR. 


•- 0.107100795 1 


STATOR LOSS COEFF. 


= 0.165054929 


ROTOR LOSS COEFF. 


•■= 0.128352524 


ROTOR LOSS THEOR. 


■- 0.418863835 


I SEN. HERD COEFF. 


_ _ _ _ 

= 4 .1 b 6 8 4 U 1 5 1^ b 1 


TH. DEG. OF REACTION = 0.951715826 i 

R C T . DEG. OF R E R C T 1 0 N = 0 . 9 8 9 8 1 5.9 5 6 : 


VELOCITY TRIANGLE 


DATA ■ 


VI = 1280.67141 


V2 = 234,80211 


VR1-- 326.98646 


VA2= 234.79915 


VU1= 1238.22418 


VU2= -1.17934 


ALPHA 1= 75.2972 


0 ALPHA 2= -0.28778 


BETA 1 = 60.6924 


3 BETA 2 = -70.61495 


POINT #8 




PRESSURE RATIO 


=•■= 2.974737893 


REF FLOW RATE 


- 1.020000000 


REF RPM 


= 17058.70933 


REF ROTOR MOMENT 


= 137,689044 


REF HP 


= 37. 2 b 7 4 9 7 9 6 


EFF. T-S 


= 0.774777057 


STATOR LOSS THEOR 


. = 0.107100795 


S T A T 0 R L 0 S S C 0 E F F 


. 0.189140155 


ROTOR LOSS COEFF. 


= 0.083159985 


ROTOR LOSS THEOR. 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ . 


= 0.414274730 


ISEN. HEAD COEFF. 


= 4.169275243 


TH. DEG. OF REACTION = 0.033929949 
ACT. DEG. OF REACT ION== 9.996859972 


VELOCITY TRIANGLE 


DATA 


VI = 1200.60325 


V2 = 219.93553 


VA1= 298.71298 


VA2= 208.69610 


VU1=- 1162. 84940 


VU2= 69.40873 



R L P H Fi 1 = 7 5 „ 5 9 3 3 4 R L P H fl 2 = 1 3 „ 3 9 6 2 4 

B £ T fl 1 = 5 9 n 5 1 8 71 BE T R 2 -79.71 2 9 8 
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TABLE III TURBINE TEST RESULTS (CONT. ) 



POINT #9 



f“‘ R E. S o 1 J R E R hi 1 1 L' 




2.519408034 


REF FLOW RRTE 


.r 


1 . 020000000 


REF RPM 


•lU 


17170,86109 


REF ROTOR MOMENT 




113.7867884 


REF HP 


=: 


3 1 . 00049275 


EFF, T-S 


rr 


0. 743362897 


STRTOR LOSS THEOR. 




0. 107100795 


STRTOR LOSS COEFF. 


:.■= 


0, 189633321 


ROTOR LOSS COEFF,, 


r: 


0. 191292857 


ROTOR LOSS THEOR. 




0.41 3 0 6 y y 4 4 


1 ci E N . HER iJ C U L F F . 




= y „ 5 6 7 6 'i! 6 1 



TH, DEG. OF RERCTIOH = O., 058183E78 
R C T . D E G » 0 F R E R C T 1 0 N = 0 „ 1 0 8 3 0 7260 



VhLOC 


ITY ’ 


■R I ANGLE 


DAT 


R 








VI = 


1 103. 


‘4 6 9 1 3 


V2 


= 


u 


+ 9.7“ 


5075 


Vfl 1 == 


259. 


39382 


VR2 


= 


ii: 


I* ’*3' C* ' 


^024 


VU1 = 


1072. 


. 54789 


VU2 




1 6 7 a y •" 


5088 


RLPHR 


1 = 


1' 6 . 4041 


r. 


ALPHA 2= 


,-i 

- 't i;« 


BETA 


•j 


57.9474 






E 1 H d - 


~ — 6 9 


POINT #10 














PRESS 


iIJRE 


RRT 1 0 


::::: 




01 


8415 


708 


k L F i" 


LOW 


RRTE 


r: 


l] 


02 


0000 


000 


REF R 


:PM 




n= 


17 


15 


8 . 62 


871 


REF R 


:OTOR 


MOMENT 






II 2 


4 6 y 8 


045 


REF HP 




r: 


cL ill 


•”1 


9 1 68 


077 


EFF. 


T-S 




rr 


0 - 


b C' 


8 1 80 


9 6 y 


STRTC 


iR LO 


SS THEOR 


- 


0 - 


10 


7100 


795 


STRTC 


iR LO 


SS COEFF, 


II *~" 


0» 


1 1 


4113 


509 


ROTOR 


: LOS 


S COEFF. 


= 


0„ 


48 


9230 


80S 


ROTOR 


; LOS 


S THEOR. 




0 » 


35 


1 6 t'l 4 


730 


I SEN. 


HERD COEFF. 








“7 I~| "7 cr 

i‘ C' i’ -J 


22027 


TH. riEG. 


0 F Rl E fl C 1 


ION 


r: 


0 . 


1147 


“i *^1 ~i •“*1 ~i 

ill ill iL_ 


RCT. 


Il L G . 


0 F R E R C 


TION=- 


•1. 


156 '9 


7E-03 


VELOC 


::ITY 


TRIANGLE 


riRTFl 








VI = 


Ci Cl "!’•* 
-• 1 


'. 09080 


l||l •-> 
1 iL. 






;45.2 


1705 


VR1 = 


2 1 0 


i„ 24041 


Vfl: 


2 = 


:i 


ir -“i ri 
■ J C' II O 


58 1 5 


VU1 = 


9 b' 


-.44140 


VLi; 


2 — 


3 


!09 a 0 


3475 



60019 
8 ’.3 f j 2 5 



FlLPHR i- 
BETR 1 



i‘ i ' n r’ U ic! o r' 

55 « 6 1 335 



RLPHR 
BETR ; 



i;;:« « .j 

-6b „ ?'} 



rdo 1 
' 6 5 9 
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t’lGURE 2 TEST RIG PIPIMG INSTALLATION 




THE FLOATING 
STATOR ASSEMBLY 



SCALt (INCHES) EOUALIY Sf-ACEO 






FIGURE 5 TURBINE C 
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FIGURE 7 TURBINE TEST RIG GEOMETRY FOR 
TURBINE CONFIGURATION C 
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FIGURE 8 ILLUSTRATION OF SMOOTHING USING 
FIRST DEGREE POLYNOMIAL 
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FIGURE 9 STATOR LOSS VS. ISENTROPIC HEAD 
COEFFICIENT, BEFORE SMOOTHING 
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2.S 3.S H.S S.S E.S 7.S B.S 

I5ENTRDPIC HEAD COEFFICIENT 






FIGURE 10 STATOR LOSS VS. ISENTROPIC HEAD 
COEFFICIENT, AFTER SMOOTHING 
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2.S 3.S H.S 5.5 E.S 7.5 B.S 

ISENTROPIC HERD CDEFFICIENT 






FIGURE 11 ROTOR LOSS VS. ISENTROPIC HEAD 
COEFFICIENT, BEFORE SMOOTHING 
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2.S 3.S H.S S.S B.S 7.S B.S 

I5ENTRDPIC HERD COEFFICIENT 



a.SB 



FIGURE 12 ROTOR LOSS VS. ISENTROPIC HEAD 
COEFFICIENT, AFTER SMOOTHING 
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FIGURE 13 PREDICTED PERFORMANCE RANGE, 

RPM=15,000, INPUT ROTOR L0SS=.251^ 
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FIGURE 14 PREDICTED PERFORMANCE RANGE, 

RPM=18,000, INPUT ROTOR L0SS=.25l4 
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FIGURE 15 PREDICTED PERFORMANCE RANGE, 

RPM=18,000, INPUT ROTOR LOSS=.l 
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FIGURE l6 PREDICTED REFERRED HORSEPOWER VS. K. 
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FIGURE 17 PREDICTED EFFICIENCY VS. 




PREDICTED EFFIECIENCY 
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FIGURE 18 PREDICTED STATOR LOSS VS. K- 
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0.70 0.72 0.7H 0.7E 0.7B 0.B0 0.B2 0.BH 0.BG 0.BB 0.30 



0.3H 



FIGURE 19 PREDICTED ROTOR LOSS VS. K. 
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EEI.eiB 



FIGURE 20 REFERRED HORSEPOWER VS. PRESSURE RATIO 
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FIGURE 21 REFERRED HORSEPOWER VS. PRESSURE RATIO 
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APPENDIX A 



TURBINE PERFORMANCE ANALYSIS 
A-1 INTRODUCTION 

This appendix gives a detailed description of the 

theory and analytical technique used to predict the 

performance of a turbine in the Turbine Test Rig. Variables 

used in this appendix are defined as they are introduced. 

Also given are the details of the computer program used 

to implement the analytical method. The variables used 

in the computer program are listed and defined in Table A-I. 

% 

A-2 ASSUMPTIONS 

The following assumptions were made* 

1. The stage was choked in the stator. 

2. In solving the continuity equation between the 
stator and rotor the solution corresponding to the lower 
static pressure was taken. 

Assumption (1) was made on the basis of past exper- 
ience with the test rig. Ref. 1 states that there is 
supersonic flow in the stator and an examination of the 
pressure distribution in the stator passages confirms this 
statement. Fig. A-1 is a plot of the pressure distri- 
bution through the stator by pressure tap number (See Fig. 6). 
It can be seen in this plot that there was a sharp pressure 
rise between tap #3 and tap #4, which indicates the 
presence of a shock in the divergent section of the passage. 
Also, the level of static pressure at the throat taps 



65 




f 



in comparison to the supply pressure, was observed to be 
independent of downstream pressure for all operating con- 
ditions . 

Assumption (2) limits the range of possible solutions 
as shown in Fig. A-2. This figure is a plot of flow 
coefficient { 0 ) vs. pressure ratio various 

loss coefficients (z). (The flow coefficient, a non-dimen- 
sional mass flux, is introduced in Section A-2.) There 
are two possible pressures for a given flow coefficient 
and given loss coefficient. Solutions were limited to 
the left side of the line connecting the points of maximum 
flow function. The procedure was chosen because it resulted 
in the "supersonic root" for the flow from the stator. 

Also, conditions corresponding to choking were sought 
for the rotor, and the locus of locally sonic conditions 
is as shown in Fig. A-3. It should be noted that because 
of the definitions of the flow coefficient and loss coef- 
ficient, the "choking" condition at which the downstreaun 
pressure has no effect upstream of the plane in question, 
does not correspond to the maximum mass flux from given 
stagnation conditions. 

The analysis is pseudo-l-dimensional. Blockage factors 
are introduced with the physical cross sectional area to 
account for non-uniform flow conditions. Loss coefficients 
are defined on the basis of kinetic energies. The analysis, 
which requires only mass flow continuity through the stage, 
is divided into three stages; the stator, the interblade 
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region and the rotor. 

While the analysis itself is general, solutions can 
only be obtained by assigning values to chosen parameters. 

In the solutions reported here,- the flow angle at the stator 
exit {<^i) , the blockage factors for the stator throat 
(k^g) and rotor exit plane and the rotor passage 

loss coefficient (Zj^) were given values. The value ofO(^ 
was fixed ato(^ =75° based on the results of the turbine 
tests and from the blading geometry. The blockage factors 
were set at k*^g=.965 and k^j^=,85. The rotor passage loss 
coefficient (zj^) was an input variable which was set at 
. 251 ^ (theoretical value) or .1 (lowest value obtained in 
test results). The "rotor loss coefficient" output from 
the program, i^^cludes all losses from the stator exit 

plane to the rotor exit plane, corresponds to the 

"rotor loss coefficient" evaluated from test rig meastire- 
mentsi consequently, 
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A-3 BASIC RELATIONS 



The sketch 
shows a general 
adiabatic process 
with entropy increase 
(losses) on a T-S 
diagram on which 
the temperature 
scale has been divided 




by the constant stagnation temperature of the process. A 
perfect gas is assumed. (0) represents the initial stag- 
nation state and ( ) represents the final state. 

The non-dimensional velocity, X, is defined as 



X = — where =\/2C It^is the "total" or "limiting” 

^to ^ 



velocity, and the loss coefficient, z, is defined as shown. 
Note that the stagnation velocity is constant throughout 
the process. 

The "flow coefficient", 0 , is defined here as 



0 = — ^ 



A(l) 



where w is the flow rate, density at the initial 

stagnation temperature and pressure, A is the flow area and 
k^ is the blockage factor. The flow coefficient defined 
in this way is a non-dimensional flow rate per unit area, or 
mass flux, referred to initial stagnation conditions. 

The flow rate is given by 
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so that Eq. A(l) can be written 




which, in terms of the non-dimensional velocity defined 
above, becomes 




The loss coefficient is defined, as shown in the above 
sketch, as 



^ = T -Tu , 

Tu-Tl, 

Rearranging Eq. A(3)» 



AD) 




/- 



XL_ 



and using the isentropic relationship between pressure and 
temperatures 




A(4) 



Using Eq. A(4), Eq. A(2) becomes a single equation 
for 0 3 iS a. function of pressure ratio analysing 

the tiirbine, generally the flow rate is known so that the 
flow coefficient can be calculated using Eq. A(l). On 
specifying a value for the loss coefficient, '( z) , the pres- 
sure ratio can be calculated from Eqf. A(2) using Eq. A(4), 
An iterative technique is used (Newton's Method) starting 
with an initial estimate of the pressure ratio. Fig. A(2) 
shows as a function of P/P^q for various values of z , 
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I 



4 



in graphical form. The solution to the left of the maxima 
is obtained by beginning the iteration with a small value 

The analysis of the turbine begins by specifying that 
the flow rate is set by the choking of the stator nozzles. 
The pressures at successive stations are then calculated 
in turn from the flow coefficient as described above. 

A -4 STATOR EXIT CONDITIONS 

The flow coefficient at the throat of the stator is the 
value of the flow coefficient corresponding to locally 
sonic conditions at the minimum area. In the present work, 
the maximum flow coefficient at a loss coefficient equal to 
.05 (defined as 0 *) was taken at the stator throat, 
throughout the calculations. 

Since 



0 * = 



W 



A(5) 



where A^ and k^^ are the area and blockage factor at the 
throat station, the flow rate being constajit through the 
machine requires that 



W = Vi* kgs - constant 



at stations ahead of the rotor, where T. is constant I this 

to 

implies that 

Ptff ^ constant 
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The value of 0 * can be obtained analytically when the 



loss coefficient is known. The maximum occurs where 

= 0, which leads to the solution (denoting values at the 

^67.) 

stator throat by an asterisk), 



( X *)2 



B 




where B = ^ 




A(6) 



A(7) 



Then 0 * is given by Eq. A(2), 

With 0 * known, continuity of flow rate requires that 
at the stator exit. 




A 

/I, Co€ cK, 



A(8) 



where subcript 1 represents the stator exit plane (station 1) 
and 0^1 is the flow angle at thestator exit. 

Knowing 0 ^^ ^i^^to found by iteration as 

described in section A-3« Solutions on the left side of 
Fig. A-2 were selected by beginning the iteration with 
close to zero. 

After obtaining thevalue for then can be 
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calculated using 



X = 



/(1-.3) 



/-I 



liL 

7 



Pu 



A(9) 



Then, found from 



Pl/P^l = (l-xj)*- 



A(10) 



and 






A(ll) 



A-5 INTERBLADE SPACE 
The interblade 
space, from Station 
(1) to Station (i) , 
is shown in the sketch. 

Conditions at (i) 
are calculated in a manner - J~J ~J ~J 
similar to those at Station (1). 

^ First, by continuity. 




irm 



0, = 0 * 

^ Pt, 






to 



A(12) 



where C{l is the flow angle at the rotor entrance. 
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The value of c<l is obtained from the condition that ang- 
ular momentum is conserved in the interblade space. Then 




A(13) 



where Rj^=mean radius at Station 1, and R^=mean radius at 
Station i. 

Then 



where N is the RPM, and the tangential velocity 
has been non-dimens ionali zed as 



Knowing , a value for P./P+^ can be found if a 

JL X X X 

value of the interblade loss coefficient (z^) is assumed. 
Then 




A(14) 



and, also 




A(15) 




A(17) 
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and 




With the absolute flow properties established at 
the rotor entrance, the relative flow conditions can be 
calculated. First, from the geometry of the velocity 
diagram, the relative flow angle ) is given by 




tan 



S JA/ oii “ LJi 

Xt CC£ 



A(19) 



where 



D. 

1 




A(20) 



is the non-dimensional rotor speed at radius R^. Then the 
non-dimensional relative velocity (X . ) is, from the vel- 

W1 

ocity diagram, given by 



^ r: 

wi “ cos 



and the temperature T^, by 



"to ^ 



A(21) 



A(22) 
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The equivalent temperature {Ref. 1), is given by 




2 



1 



A(23) 



and therefore 







T. / P 



to 



A(24) 



where 

Tgi = equivalent temperature into rotor (See Fig. A-7) 
Pgi = equivalent pressiare into the rotor (See Fig. A-?) 

a-6 rotor conditions 

Using the values calculated thus far it is now possible 
to calculate the conditions in the rotor. Continuity 
requires that 



where a subscript R denotes the exit plane in the rotor 
blading. 

It should now be noted that the process in the rotor 
frame from the equivalent conditions at (i) to the exit 
of the rotor is entirely similar to the basic process 
described in Section A-3* However, the equivalent temper- 
ature takes the place of stagnation temperature and 
relative velocity replaces velocity in the given equations. 
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The pressure ratio which satisfies the flow coefficient 



given by Eq. A(25). for an assumed rotor passage loss 

This is obtained as described 
in Section A-5* The corresponding non-dimensional velocity 
is now 



coefficient (zj^) is 




A(26) 



where 



W, 



w2 






A(27) 



Hence 



'T„. 

Y = Y / - El 

^w2 w2yT!^ 



A(28) 



From the velocity diagram, 



4 = <^^2 - ^ W 2 * (=' w 2 



A(29) 



where = 71° (from the geometry of the rotor), and 
— _ Ro 

"2 = «i-Rf • 

The temperature at the rotor exit is given by 



- 2 - - 



, JLL ‘ 


/ \ 


a-H) 







to 



^t2 ^2 2 

so that ~T — = T^ 
to to 



A(30) 

A(31) 
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Next the loss coefficient that includes all of the 
losses in the stage from Station 1 to Station 2 can be 
calculated. This coefficient includes all of the losses in 
the inter-blade region as well as those in the rotor 
passage. It is called z ^_2 here but, in fact corresponds 
to the measured rotor loss coefficient in the TTR since all 
losses aft of the stator in the TTR are included in the 
"rotor loss coefficient." 



, = 1 . 

^1-2 ^ T„-T, 

^ ns 



H TL. 




n, Ttc 
_2k _ZL 

= 1 - 21 . 

..TV.,.. 

Tct to 

A-7 STAGE PERFORMANCE PARAMETERS 



A(32) 



In the notation used for the turbine test rig in 
Ref. 1 and Ref. 2, the following equations determine the 
turbine performance parameters from quantities calculated 
in the above steps* 



22vT 



w 





^^is 





A(33) 

A(34) 




T-S 



ATis 



A(35) 
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A(36) 



H.P. 

M = 

cT = - 
e = 

w* = 
H.PJ 
M* = 
N* = 



= w (.2402) 



A(37) 



H.P. (550) 1^ 

frN 





w/r 

"T~ 



A(38) 

A(39) 

A(40) 

A(4l) 



= H.P. 




A(42) 

A(43) 









A(44) 
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a-8 the computer program 

A-8.1 Overview 

There are nine degrees of freedom in the com- 
puter program* 

1. CP^f , the flow angle out of the stator 

2. 1 ^ 13 * "the stator exit blockage factor 

3. interstage blockage factor 

4. "the rotor passage blockage factor 

5* Zg, the stator loss coefficient 

6. z^, the interstage loss coefficient 

7. Zj^f the rotor passage loss coefficient 

8. N, turbine RPM 

9. "total temperature upstream of the stator 

c<, is input on the basis of turbine test results. 

At high pressure ratios the value of was measured to 

be 75°-76°. Therefore, the value of o<, was set and held 

at 75° • The values of RPM and T^^ are tho® conditions 
for which a solution is being sought. 

All the other parameters, the rotor passage loss 
coefficient (z^^) and blockage factor were given 

chosen values. This choice was made because the incidence 
losses were included in the interblade calculation. It 
Was thought that the parameters describing the flow inside 
the rotor passages could be held fixed as the speed 
was allowed to vary. 

Consequently, there are four degrees of freedom with 
which the program works. The goal for the program is to 
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find those solutions which result in flow coefficients in 
the rotor that are less than, or equal to, the "choking" 
value. This is done in the program by iteration on the 
remaining four degrees of freedom. Fig. A-5 is a block 
diagram showing schematically the iteration process. The 
order of iteration is (from the most frequent to the least 
frequent), , z^, and The loss coefficients for 

both the stator and interblade always begin at zero and they 
are incremented as the program seeks solutions. A graphical 
representation of a typical iteration is shown in Fig. A-4. 

It can be seen that the program can find a solution for 
the condition where the loss coefficient is equal to zero 
(denoted as (1) on Fig. A-4). On the next iteration the 
program will set the loss coefficient equal to .1 and it 
can be seen that a solution still exists at point (2). On 
the next iteration, however, it can be seen that a solution 
does not exist for a loss coefficient equal to .2. In this 
case the program will go back and add a smaller increment 
to the previous loss coefficient (.1). 

Each time an added increment results in a loss coef- 
ficient greater than the maximum allowable the program 
subtracts that increment and adds a smaller one. In this 
way the program eventually finds the maximum loss coef- 
ficient which will yield a solution. At this point, the given 
flow coefficient is at the maximum point for the calculated 
loss coefficient line. At this point "choking", as defined 
here, has occurred. 
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It should be pointed out that this "choking” con- 
dition (the point of maximum flow coefficient for any 
given loss coefficient) is not the point corresponding 
to locally sonic conditions. In Fig. A-3 the locus of 
points of maximum flow coefficient and the locus of points 
of local sonic flow are both shown. The point of locally 
sonic flow always falls to the left of the "choking" point 
(as defined in this work), and the flow coefficient for 
sonic flow is always less than the flow coefficient at 
the defined choking point. The reason is that the point 
of maximum mass flux (maximum flow coefficient) is for 
given upstream stagnation conditions (see next section). 
Because of the way the flow coefficient is defined, the 
maximum for a given loss coefficient will not correspond 
to local sonic conditions arxcept where the loss coefficient 
is equal to zero. 

The difference between the "choking" flow coefficient 
corresponding to the maximum and the flow coefficient at 
the local sonic condition is about 10^ at a loss coefficient 
of .3 (the highest rotor loss coefficient predicted by 
the prograjn) . The difference is about at a loss coef- 
ficient equal to .25 (the highest stator loss predicted), 
and insignificant (less than \fo) at a loss coefficient of 
about . 15 . 
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A-8.2 Description of the Program 

The following discussion refers to Fig. A-6, 
which is a detailed flow chart of the program, to 
Table A-I , which is a list of program variables and their 
definitions, and to the program listing given in Section A-10. 

The first step in the program is to input the variables. 
o(f is input in line 280 (it is also converted to radians 
in that line), is input in line 286 and k^^^ is input 

in line 284. All other variables are input as requested 
by the program (lines 20-110). Note that and are 

input also. They are not involved in the calculations and 
were originally input as the conditions for which a 
particular solution was sought. 

The next step is to go to a subroutine to calculate 
the value of (line 1?0). j6* is defined as the value 

of the flow coefficient at the throat of the stator which 
sets the values of the flow coefficients at all downstream 
stations. is generated in the subroutine by the method 

covered in Section A-4 assuming a loss coefficient of .05 
to the throat of the stator (line l 60 ). 

Next the program calculates the conditions in the 
stator exit plane. This is done by first calculating the 
value for (line 320). As explained earlier, the value 
for determines the upper limit of the stator loss 
coefficient. 

Using the method described in Section A-3, the program 
calculates the value for Note that if the stator 
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loss coefficient is large enough, so that a solution does 
not exist, this will he detected in the decision, "P^/P^^<0 
or>l" (lines 350 and 370). If "yes" then the stator loss 
(z_) is too large and the following steps are takeni 

O 

1. The last increment added to the stator loss is 
subtracted (line 420). (Note, the initial stator loss 
increment is .1.) 

2. The existing stator loss increment is multiplied 
by some factor less than 1 (line 430), 

3* The new, smaller, increment is then added back 
to the stator loss coefficient (line 1570). 

If stator loss increment is sufficiently small 
(line 1575) then the upper limit of the loss coefficient 
has been reached as discussed earlier. At this point no 
further solutions are attainable and k^^^ is incremented 
(line 2290). Note that the upper limit of stator loss 
will not normally be attained until several iterations 
have been made through the entire prograin. 

Next the interblade calculations are made. The sol- 
ution for found for the stator exit plane is used 

to make the calculations for 0 ^ (line 770)* Again, 0 ^ 
determines the upper limit on the interblade loss. In 
the interblade region, however, there is another criterion 
that must be met. The value for generated within the 
program must be equal to (^^ =69° and is fixed by the 
geometry). This criterion is met by iterating the inter- 
blade loss coefficient as follows: On every pass through 
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this section the interblade loss is initially set to zero 
and the interblade loss increment is initially set to .1 
(lines 720, 7^0). The value for is calculated using 

the same procedure as described earlier for Again, 

the value of is checked on each iteration to ensure 

that it falls between 0 and 1 (line 780). If it does not 
then the interblade loss is reduced in a method analagous 
to that covered for the stator loss. 

After generating a value for Pj^/P-^j^ "the program cal- 
culates (line 1208). When the value for is generated 
it is compared to and three results are possible. 

If is equal to ^^+.05 then converganne has occurred and 
the prograjn will continue. If is less than (line 1250) 
then the value for interblade loss coefficient is increased 
by the current value of interblade increment and the program 
goes back to the start of the interblade calculations and 
begins anew. If is greater than (line 1240) the 
value of interblade loss is decreased in the same way 
as was described for the stator loss. 

Two results are possible in the interblade regiom 
Either the program finds the value of the interblade loss 
that gives flow angle convergence ( ^; “ ^ inter- 

blade loss increment becomes sufficiently small (10“^) to 
trigger termination. 

Note that the interblade loss increment is allowed to 
become very small (10“^) (line 850). This is because the 
above mentioned angle convergence requirement is very 
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sensitive. Note, also, that when the interblade loss 
increment falls to lO”^ and the iteration process in the 
interblade region is terminated, the stator loss coefficient 
is incremented (line 850 ). Stator loss is incremented 
because a higher value of stator loss may result in con- 
ditions that allow convergence in the interblade region. 

If not, the stator loss calculations eventually reach the 
termination criterion mentioned earlier. 

Once angle convergence has been attained in the inter- 
blade region then the value of the rotor flow coefficient 
can be calculated (line 1470). If the program gets to this 
point for any given set of conditions then a solution 
exists and the only determination to be made is whether 
the resulting rotor flow coefficient is greater than the 
••choking” value (line 1530). If ”yes" then it is dis- 
allowed, since a flow coefficient greater than the "choking" 
value cannot occur physically. If the calculated coef- 
ficient is less than, or equal to the "choking" value then 
the results are printed, the stator loss is incremented, and 
the calculations start again at the stator exit plane. 

The choking value of the flow coefficient for the 
rotor (0*^) Is determined in exactly the same way as the 
choking value for the stator (0*)t with the assumed value 
of rotor loss used as the loss coefficient. (At present 
the value for 0^ is calculated separately and inserted in 
line 121. It is suggested that a "GOSUB 2l40" be inserted 
at this point. Then the input value for rotor loss 
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coefficient can be changed easily and the “choking" value 
for the rotor will be calculated in the program.) If 
the resulting flow coefficient is greater than "choking", 
the program increments stator loss coefficient and returns 
to the stator plane calculations to begin anew. 

Therefore, in all cases the program will seek the 
highest value of stator loss coefficient for which solutions 
exist. At some point, the rotor flow coefficient normally 
attains a value less than the "choking" value and all 
solutions for higher stator, loss coefficients will be 
acceptable solutions. Or, if the "choking" condition is 
never reached, the program will attain the highest possible 
stator loss in attempting to achieve solutions less than 
"choking". 

When the highest value of the stator loss is reached, 

there are two possible avenues in the program. If the 

program was generating solutions less than "choking" at 

the time of termination in the stator calculations, then 

kbi will be incremented, stator loss will be set back to 

zero, and the process begun anew (line 2280). If the 

solutions being generated at the time of termination were 

/ 

greater than "choking", then the upper value of k^^ has 
been reached and any increase in k^j^ will yield no further 
solutions. In this case k^^ is incremented, k^^ is set 
to the pre- determined lower limit, and stator loss is 
reset to zero (line 2330 ). 

The program continues in the above manner until the 
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pre -determined upper limit of is reached (line 2335 ) i 
at which point the program is terminated. 

Since the 21-MX computer has only 6 digits of accuracy, 
round-off error has been a problem. Where Newton's Method 
is used to solve for pressure ratio the computer lacks 
sufficient accuracy to converge to the desired accuracy 
under certain conditions. When this happens, the program 
will stay inside the iteration and never converge. This 
problem has been solved by putting a counter inside each 
iteration loop (lines 5^0 and 980). More than 60 iterations 
is treated as a non-convergent condition, and the approp- 
riate loss coefficient is reduced in the manner covered 
previously. 

A-9 OPERATING PROCEDURE 

The procedures for operating the Hewlett-Packard 
21-MX computer will not be covered in this paper since 
the applicable manuals are available at the laboratory. 

The computer must be on, with the RTE-B operating system 
••READY". 

1. Load the paper tape program labeled "TTR 11". 

2. Edit into the program the minimum values of k^^ 
(line 284) and k^^ line (286). Edit the desired value of 
cK, into line 280. If it is desired to get a print-out 

of every solution then either remove lines 1724 and 
1726 or place ••REM" in front of them. If these two lines 
are left in the program then a solution will be printed 
out only if the calculated pressure ratio (P2^^^to^ less 
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than, or greater than, the previously calculated lowest or 
highest pressure ratio respectively, for the value of 
then being used. 

3. Type "RUN”, "RETURN" on the keyboard and the 
program will begin execution by asking for inputs? 

"INPUT STATOR LOSS COEFFICIENT” The normal input is 
zero. However, if some specific case is required then any 
desired loss coefficient can be input. 

"INPUT PTO” Input the upstream tetal pressure. 

"INPUT P2” Input the hood pressure. 

"INPUT RPM" Input the RPM desired. 

"INPUT TTO" Input the total temperature into the 
stage . 

This is all that is required to operate the program. 
Depending on the range of k^^ requested the program will 
take from 30 minutes to 36 hours to run on the 21-MX. 

It is important to realize that the above procedure 
will produce solutions that give values of the rotor 
flow coefficient less than "choking". If the turbine 
performance at the "choking" chndition is desired then 
after the computer run is complete a further step must be 
taken. It is necessary to take the values of k^^ and k^^ 
for which a solution is desired and force the choking 
condition. This is done as follows? 

1. Scan the output results at the desired k^^ and 
k^^ and locate the point that has a value of rotor flow 
coefficient closest to the choking value and a non-zero 
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stator loss coefficient. 



2. Re-start the program, this time putting in a value 
of the stator loss coefficient slightly less than the value 
printed on the output. 

As stator loss is manually decreased, the rotor flow 
coefficient will increase towards the choking value. When 
the rotor flow coefficient is within the limits specified 
in line 2370 the words "CONVERGENT CONDITIONS" will appear 
on the teletype output preceding the printed results. 

If the stator loss coefficient is decreased too much, 
the rotor flow coefficient will increase to a value above 
choking. In this case no ouput will be printed on the tele- 
type. However, the flow coefficient will appear on the 
video display. In this case increase the stator loss 
slightly. 

It has been found in this work that the extreme values 
of referred horsepower, stator loss, rotor loss, and 
efficiency occur at the lowest and highest values of k^^^ 
for each value of k^^. With this in mind it is possible to 
find the range of predicted values at choking by forcing 
to choking only the lowest value of k^^ (with non-zero 
stator loss) and the highest value of k^^ for each k^^. 
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2A5 REM INPUT AHt.A OK STATOR THROAT 
230 LET A0=. 020625 

255 REM INPUT bLOCKAGE FACTOR h OR STATOR 
260 LtT Kl=.965 

265 REM INPUT AnnA OF INTERSTAGE RcoION 
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1300 L£.T ZS=Z9*r.8 
1305 IF Z3<0 TH£iN 2270 
1308 IF Z9<. 000001 Th£N 860 



1320 PRINT# r*AT ANGLE CONVERGENCE KOLLOWING CONDITIONS EXISTED 
1330 PRINT# r*STATOR LOSS COeKFICIENT="eO 
1340 PRINT# 1 "INTeRSThG£ LOSS COeFFI CIENT = "E 3 
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2330 LET K6=K6+.01 



2335 IF (K6>.«2) TH£N 2620 

2336 LET Ui8 = 0 

2337 LET 

2340 LET K4:.4 
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TABLE A-I 



AO 

A1 



DEFINITION OF VARIABLES 

ik*) Area at the Stator Throat 
(A^) Area at the Stator Exit Plane 



A8 


Decision Variable 




A9 


Decision Variable 




Cl 


Decision Variable 




D 


(^>’Ptc/Pref 




D1 






D2 


H-kf 




D3 






D4 






E 


Decision Variable 




EO 


(7t_s^ Total-to-Static 


Efficiency 


FO 


Flow Coefficient 


at Stator Throat 


FI 


(0^) Flow Coefficient 


at Stator Exit Plane 


F2 


(0^) Flow Coefficient 


Through Rotor Passage 


F3 


(0^) Flow Coefficient 


Through Interblade Area 


F4 


Input Value of 0^ Based on Assumed Rotor Loss 


F8 


Decision Variable 




F9 


Decision V-^riable 




GO 


(Gamma - 1)/Gamma 




HO 


Horsepower 




HI 


Referred Horsepower 




H8 


(5*) Rho 




H9 


Reynolds Number 
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J 



{rr ) Pi 



K1 Blockage Factor at Stator Throat 

K4 k^^ Interhlade Blockage Factor 

K5 k^j^ Blockage Factor in the Rotor 

K6 k^j^ Blockage Factor at the Stator Exit Plane 

Ml Rotor Torque 

M2 Referred Rotor Torque 

N RPM 

NO Referred RPM 

01 (^i) Flow Angle at Stator Exit Plane 

03 i^i) Flow Angle at Rotor Entrance Plane 

04 (Ak^^) Increment for k^^ 

'’to 

P' Pl/Pto 

P2 P/Pto 

P3 Pl/Ptl 

P** Ptl/Pto 

P5 Pi/Ptl 

P6 Pi/Pto 

P7 PEi/Pto 

P8 Decision Variable 

P9 P 2 

F (jto) Total Density at Stagnation Conditions 

R7 Relaxation Parameter 

R8 Relaxation Pair ame ter 

R9 Predicted P2/^to 
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S7-S9 

TO 

T1 

T2 

T3 

T4 

T5 

T6 

U1 

U2 

U3 

VO 

WO 

W1 

W9 

xo 

XI 

X2 

X3 

X4 

X5 

X6 

X8 

Y2 

Y6 

Y7 

Y8-Y9 



Used in GOSUB 7000 to Find P2/^Ei 

V^to 

Dimensionless Velocity (See Sppendix A) 
Dimensionless Velocity (See Appendix A) 
Dimensionless Velocity (See Appendix A) 
V^^ Total Velocity 
w Flow Rate 

Referred Flow Rate 
Decision Variable 

X^^ Non-Dipiensional Relative Velocity 
X^ Non-Dimensional Velocity 
X 2 Non-Dimensional Velocity 
X^ Non-Dimensional Velocity 
Xui Non-Dimensional Velocity 
X . Non-Dimensional Relative Velocity 

W1 

X ^2 Non-Dimensional Relative Velocity 

Decision Variable 

Non-Dimensional Velocity 

Increment for 

Increment for ^bi 

Counters 
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ZO Zg Stator Loss Coefficient 

Z3 Interblade Loss Coefficient 

Z7 ^1-2 Predicted Rotor Loss Coefficient 
Z8 Increment for Stator Loss Coefficient 

Z9 Increment for Interblade Loss Coefficient 
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GCI.CIB 



FIGURE A-1 PRESSURE DISTRIBUTION THROUGH THE STATOR 
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B.eiei i.Bia 2.m a. cm h.bci s.bb s.em 

PRESSURE TRP ♦ 



FIGURE A-2 ILLUSTRATION OF SOLUTION AREA 
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P I /PT0 



FIGURE A-3 MAXIMUM 0 AND SONIC LINE 
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FIG. A-5 PROGRAM SCHEMATIC 
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FIG. A-6 flow chart 
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FIGURE &-?. THERMODYNAMIC PROCESS OF FLUID 

IN AN AXIAL TURBINE STAGE 






APPENDIX B 



TURBINE TEST RIG (TTR) DATA REDUCTION AND PROCESSING 
B-1 INTRODUCTION 

This appendix describes only those changes that have 
been made to the data reduction procedure given in Ref. 2. 
The prograjn numbers and the functions of those programs 
given in Ref. 2 have not changed. The channel and port 
assignments for data collection and storage have not 
changed. Variables also have not changed, although 
additional variables were defined in the coxurse of mod- 
ifying the data reduction process. Those variables that 
now exist in addition to those given in Ref. 2 are given 
in Table B-I. 

Modifications made during the course of this work were* 

1. Raw data can now be read directly from mass memory. 
Previously it was necessary to read the paper tape for a 
given point in order to reduce it. 

2. Provisions have been written into the programs 
for smoothing the reduced data. 

3. All of the parameters used in calculating the 
theoretical loss coefficients that were previously 
entered from charts are now in polynomial form in the 
program. 

4. Storage of raw data is now a separate program. 
Previously it was necessary to run the reduction sequence 
to store data. 
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B-2 DESCRIPTION OF PROGRAMS 

The data reduction is divided into ten separate 
programs with an additional program to store the raw data. 
Fig. B-1 shows the contents of each program and also shows 
the sequence in which the programs are chained. Following 
is a description of the reduction sequence referred to 
Fig. B-1 which will be followed by a description of how to 
run the reduction program. 

TTR in Fig. B-1 is used only to store the raw data. 
TTR does not chain to any other program. 

The data reduction process begins in TTRl. In order 
to smooth the reduced data as discussed in Section III it 

is necessary to have n values for (Pi“Phub^^^^tip*’^hub^ * 
Therefore, the first step in the data reduction process is 
to go from TTRl to the point in TTR2 where the values of 
Pi , and have been calculated. On completion, 

the program returns to TTRl and repeats the process n 
times until the n values of (Pi~Phub^'^^^tip'’^hub^ have 
been calculated. The program then proceeds to TTR9 where 
a polynomial curve fit for the n values of 

is generated. Then the prograjn 
returns the polynomial coefficients to TTRIB where the 
reduction process begins at the first point in the run, 
However, when the point in TTR2 is reached where is 

calculated, it is calculated using the polynomial as 
described in Section III. 

The rest of the reduction process in Fig. B-1 is the 
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same as described in Ref. 2 with the exception that no 
further keyboard inputs are needed until the raw and 
reduced data is tabulated. 

B-3 RAW DATA STORAGE 

The procedure for storing raw data on the mass memory 
is outlined in Ref. 2, Appendix B. Sec. B-3t lines 10t19» 

Note that line 10 should read, ”GET *TTR"'. Although it is 
now a separate program the procedures for storing data 
have not changed. 

B-4 PROCEDURE FOR DATA REDUCTION PROGRAM 

1. Key in GET "TTRIB" 

2. Press "RUN EXECUTE". 

3. The following check list will be printed on the 
HP 9830 printer! 

"PRIOR TO RUNNING THIS SEQUENCE ENSURE FOLLOWING! 

"THAT TTR2 LINE 590 HAS PROPER FACTOR IN IT; 

1.02 FOR HOODED 
1.01 FOR UNHOODED 

"IF IT IS NECESSARY TO INPUT ALPHA 1 THEN CHANGE TTR2 
"LINE 960 TO 'INPUT A3' 

"IF BLOCKAGE FACTOR OTHER THAN 1, CHANGE TTR2-1060 
"TO 'INPUT X7', AND PUT SEMICOLON AFTER TTR2-1050 
"ENSURE THAT TTRl, TTRIB HAVE THE PROPER FILENAME IN 44l 
"IF IT ISN'T DESIRED TO STORE REDUCED DATA CHANGE 
TTR 6 - 58 O TO 'G1=0' 

"ENSURE TTR 6 - 6 IO, 630 HAVE PROPER FILENAME FOR REDUCED 



DATA 



114 



"ENSURE TTR7-220 HAS PROPER FILENAME FOR RAW DATA 

"ENSURE TTR8-220 HAS PROPER FILENAME FOR REDUCED DATA 

"IF IT IS DESIRED TO RUN A SINGLE POINT, WITHOUT 
SMOOTHING, INSERT TTR2-1127 'GOTO 1200'. OTHERWISE 
OMIT TTR2-1127 

The above checklist will prepare all ten programs. The 
last item is a provision for elimination of smoothing if 
it is desired to run a single point. If smoothing is not 
desired then each point must be run separately with 
TTR2-1127 inserted, 

4. Press CONTINUE EXECUTE 

5. The display will read LOWEST, HIGHEST RECORD # 

THIS RUN. Input the lowest record number and the 

highest record number from which it is desired to read • 
raw data. Note, the record number on which the reduced 
data will be stored will be the same number as the raw 
data record from which the raw data was read. Make sure 
the filename designated for the storage of the reduced 
data can accept data on those record numbers without 
writing over previous information. 

6. Next the display will read PRINT OUT RESULTS? 

YES=1, N0=0, If a 1 is selected here then at the conclusion 
of data reduction, for all the points in the run, the 
program will begin to print all results. If 0 is selected 
then all reduced data will be stored on the mass memory. 

7. It will take about 15 minutes for the program to 
read the raw data and make the calculations necessary for 
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s 



<3 



smoothing the data. When this is oompleted the program will 
start reducing each point and printing the results. 

8. At the completion of data reduction for all 
points the option to tabulate the data is available. 

ENTER RECORD #’S: LOWEST, HIGHEST will appear on the display. 
If it is desired to tabulate the raw data then enter the same 
record numbers that were entered to initiate the program. 

The raw data will then be tabulated. 

9. Next the display will read ENTER LOWEST, HIGHEST 
RECORD NUMBER. If it is desired to tabulate the reduced 
data then enter the appropriate record numbers. The 
reduced data will be tabulated. 
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TABLE B-I 



VARIABLES ADDED TO DATA REDUCTION PROGRAM 



A2 


Lowest Record Number 




A5 


Counter 




A6 


Decision Variable 




A? 


Highest Record Number 




B 


(Array) Polynomial Coefficients From 


TTR9 


F8 


Decision Variable 




Q 


(Array) Values of Isentropic Head Coefficient 


T4 


Polynomial Approximation of Reynolds 


Number 


Z 


Decision Variable 
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1 . ) RAW DATA STORAGE 

2. ) RAW DATA PRINTED 




FIGURE B-1 DATA REDUCTION SCHEMATIC 
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